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ABSTRACT 

Detailed Monte Carlo Inversion analysis of the spectral lines from three Lyman limit systems 
(LLS) [A(Hi) > 1.0 xlO 17 cm" 2 ] and nine lower AT(Hl) systems [2 x 10 14 cm" 2 < iV(Hl) <2x 
10 16 cm" 2 ] observed in the VLT/UVES spectra of Q0347-3819 (in the range 2.21 < z < 3.14) 
and of APM BR J0307-4945 (at z = 4.21 and 4.81) is presented. Combined with the results from 
a previous work, the analyzed LLSs show that they are a heterogeneous population originating 
in different environments. A functional dependence of the line-of-sight velocity dispersion er v on 
the absorber size L is confirmed: the majority of the analyzed systems follow the scaling relation 
cr v ~ (Ah£)°- 3 (with Ah being the total gas column density). This means that most absorbers 
may be related to virialized systems like galaxies or their halos. Previously noted enhancement 
of the metal content in small size systems is also confirmed: metallicities of Z ~ (1/3 — 1/2) Zq 
are found in systems with L < 0.4 kpc, whereas we observe much lower metal abundances 
in systems with larger linear sizes. For the first time in LLSs, a pronounced [a-element/iron- 
peak] enrichment is revealed: the absorber at z a b s = 2.21 shows [O/Fe] = 0.65 ±0.11, [Si/Fe] = 
0.51 ±0.11, and [Mg/Fe] = 0.38 ±0.11. Several absorption systems exhibit characteristics which 
are very similar to that observed in high- velocity clouds in the Milky Way and may be considered 
as high-redshift counterparts of Galactic HVCs. 

Subject headings: cosmology: observations — line: identification — quasars: absorption lines — quasars: 
abundances — quasars: individual (Q0347-3819, APM BR J0307-4945) 
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1. Introduction 

With the present paper we continue to study 
the chemical composition and the kinematic char- 
acteristics of quasar absorption systems using a 
new computational procedure, - the Monte Carlo 
Inversion algorithm (MCI), - developed earlier in 
a series of papers [see Levshakov, Agafonova & 
Kegel (2000); hereafter LAK]. The MCI technique 
allows us to recover self-consistcntly the physi- 
cal parameters of the intervening gas cloud (such 
as the average gas number density no, the col- 
land 
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umn densities for different species JV a , the kinetic 
temperature Tki n , the metal abundances Z a , and 
the linear size L), the statistical characteristics of 
the underlying hydrodynamical fields (such as the 
line-of-sight velocity dispersion ov, and the den- 
sity dispersion a y ), and the line of sight density 
nn(x) and velocity v(x) distributions (here x is 
the dimensionless coordinate in units of L). Hav- 
ing this comprehensive information we are able 
to classify the absorbers more reliably and hence 
to obtain important clues concerning the physical 
conditions in intervening galaxies, galactic halos 
and large scale structure objects at high redshifts. 
Besides, it will also be possible to constrain the 
existing theories of the origin of galaxy formation 
since the observed statistics of the damped Lya 
(DLA) and Lyman limit (LLS) systems is believed 
to be a strong test of different cosmological mod- 
els (e.g. Gardner et al. 2001; Prochaska & Wolfe 
2001). 

In the first part of our study (Levshakov et al. 
2002a, hereafter Paper I) we reported results on 
the absorption systems at z a ^ s = 1.87, 1.92 and 
1.94 toward the HDF- South quasar J2233-606. 
These systems exhibit many metal lines with quite 
complex structures. It was found that all profiles 
can be well described with an assumption of a ho- 
mogeneous metallicity and a unique photoioniz- 
ing background. According to the estimated sizes, 
velocity dispersions and metal contents the ab- 
sorbers at z a bs= 1-92 and 1.87 were related to 
the galactic halos whereas the system at z a b s = 
1.94 was formed, more likely, in an irregular star- 
forming galaxy. It was also found, that the linear 
size and the line-of-sight velocity dispersion for all 
three absorbers obey a scaling relation of the same 
kind that can be expected for virialized systems. 

The present paper deals with absorbers ob- 
served in the spectra of Q0347-3819 (z cm = 3.23) 
and APM BR J0307-4945 (z cm = 4.75, see § 2.1). 
Both spectra include several dozens of systems 
containing metals, but most of them are weak and 
severely blended and hence do not allow to esti- 
mate the underlying physical parameters with a 
reasonable accuracy. After preliminary analysis 
only 12 systems were chosen for the inversion with 
the MCI and their properties are described below. 

The structure of the paper is as follows. § 2 
describes the data sets. In § 3 our model assump- 
tions and basic equations are specified. The esti- 



mated parameters for individual systems are given 
in § 4. The implication of the obtained results 
to the theories of LLS origin are discussed in § 5 
and our conclusions are reported in § 6. Appendix 
contains a table with typical parameters of differ- 
ent absorbers which are referred to in the present 
study. 

2. Observations 

The spectroscopic observations of Q0347-3819 
and APM BR J0307-4945 obtained with the UV- 
Visual Echo-lie Spectrograph UVES (D'Odorico et 
al. 2000) on the VLT/Kueyen 8.2 m telescope 
are described in detail by D'Odorico, Dessauges- 
Zavadsky & Molaro (2001) and by Dessauges- 
Zavadsky et al. (2001), respectively. Both spectra 
were observed with the spectral resolution FWHM 
~ 7 km s _1 . For the analysis of metal systems 
from the Q0347-3819 spectrum with lines in the 
range 4880 - 6730 A which was not covered by 
the VLT observations, we used a portion of the 
Q0347-3819 spectrum obtained with the High- 
Resolution Echelle Spectrograph HIRES (Vogt et 
al. 1994) on the 10 m W. M. Keck I telescope 
(Prochaska & Wolfe 1999). The spectral reso- 
lution in this case was about 8 km s _1 . The 
VLT/UVES data are now available for public use 
in the VLT data archive. 

The majority of the metal systems in the spec- 
trum of Q0347-3819 were identified in Levshakov 
et al. (2002b), whereas the z a b s = 4.21 system to- 
ward APM BR J0307-4945 was distinguished by 
Dcssauges-Zavadsky et al. 7 (2001) as consisting 
of two sub-systems : one at z a b s = 4.211 and the 
other at z a b s = 4.218. A new system at z a b s =4.81 
is analyzed here for the first time. 

2.1. Emission redshift of APM BR J0307 
4945 

The emission redshift of this distant quasar 
z cm = 4.728 ± 0.015 was previously measured by 
Pcroux et al. (2001) from the Si IV+O IV] A1400.0 
and C IV A 1549.1 lines observed in the ~ 5 A 
resolution spectrum obtained with the 4 m Cerro 
Tololo Inter- American Observatory telescope. 

7 Data are listed in Table 5 which is only available in elec- 
tronic form at the CDS via anonymous ftp to cdsarc.u- 
strasbg.fr (130.79.128.5). 
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In our VLT/UVES spectrum of this quasar a 
few additional lines can be identified which are 
useful for the redshift measurements. The most 
important of them is the weak O I A1304 line. 
From earlier studies (see, e.g., Tytler & Fan 1992 
and references cited therein) it is known that 
'low-ionization lines' such as O I A1304 are sys- 
tematically redshifted and narrower than 'high- 
ionization' lines such as C IV, Si IV, and Lya. 

In Fig. 1 we compare the O I profile with those 
of C IV and of a wide blend of the Lya+N V+Si II 
lines. All these lines are shown in the same ve- 
locity scale which is defined by the O I A1304 
center corresponding to z cm = 4.7525. This O I 
line is redshifted with respect to the z em value de- 
duced by Peroux et al. from the measurements of 
the Si IV+O IV] and C IV profiles. Because the 
Lya emission line is blended with other emission 
lines and its blue wing is distorted by numerous 
narrow absorption lines, profile comparison can- 
not be very accurate in this case. Nevertheless, 
we found a smooth fit to the C IV A1549 profile 
(which is unblended and shows significant asym- 
metry) and used this synthetic profile for compari- 
son with other lines (altering the amplitude of the 
synthetic profile to match the line profile while 
keeping its center unchanged). 

Fig. 1 shows that this simplified procedure in- 
deed allows us to achieve a fairly good concordance 
between the C IV line and the Lya+N V+Si II 
blend. This could indicate that the redshift of the 
quasar is higher than that measured by Peroux et 
al., being actually z cm ~ 4.753. 

3. Model assumptions and the MCI pro- 
cedure 

The complete description of the MCI code is 
given in LAK and its most updated version - in 
Paper I. Since this technique is relatively new, we 
briefly outline it here and stress its difference from 
the Voigt profile fitting (VPF) procedure com- 
monly used for the QSO absorption line analysis. 

The VPF deconvolution is based on the as- 
sumption that the observed complex line pro- 
files are caused by several separate clouds ran- 
domly distributed along the line of sight. In every 
cloud, gas is characterized by the constant density 
and normally distributed velocities (the b param- 
eter usually estimated in the VPF procedure just 



stands for the dispersion of the velocity distribu- 
tion). Because of the constant gas density, the 
ionizing structure inside the cloud is described by 
a single ionization parameter U which can be esti- 
mated from the measured column densities of lines 
of different ions with help of some photoionization 
code if the spectrum of the background ionizing 
radiation is given. 

However, numerous cosmological hydrodynam- 
ical calculations performed in the previous decade 
have shown that the QSO absorption lines arise 
more likely in the smoothly fluctuating intergalac- 
tic medium in a network of sheets, filaments, and 
halos (e.g., Cen et al. 1994; Miralda-Escude et 
al. 1996; Theuns et al. 1998). This model finds 
its support also in modern high resolution spectro- 
scopic observations: the increasing spectral resolu- 
tion reveals progressively more and more complex 
profiles. A very important characteristic of the 
continuous absorbing medium is that the contri- 
bution to any point within the line profile comes 
not only from a single separate area (a 'cloud') 
but from all volume elements ('clouds') distributed 
along the sightlinc within the absorbing region and 
having the same radial velocity (see, for details, 
§ 2.2 in LAK). 

If the absorption systems are indeed formed 
in the fluctuating IGM, then the above described 
VPF procedure which interprets each absorption 
feature in the line profile as caused by one distin- 
guished cloud is not, in general, physically justi- 
fied. In some systems, this approach can produce 
rather erratic results like extremely varying metal- 
licities between subcomponents, negative kinetic 
temperatures, exotic UV background spectra etc. 
(see examples in Levshakov et al. 1999; LAK; Pa- 
per I). 

The MCI procedure is based on the assumption 
that all lines observed in a metal system arise in 
a continuous absorbing gas slab of a thickness L 
with a fluctuating gas density and a random veloc- 
ity field. We also assume that within the absorber 
the metal abundances are constant, the gas is op- 
tically thin for the ionizing UV radiation, and the 
gas is in the thermodynamic and ionization equi- 
librium. The last assumption means that the frac- 
tional ionizations of different ions are determined 
exclusively by the gas density and vary from point 
to point along the sightline. These fractional ion- 
ization variations are just the cause of the observed 
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diversity of profile shapes between ions of low- and 
high-ionization stages. 

Whereas most of the above mentioned assump- 
tions are quite natural, that of the constant metal- 
licity over the entire absorbing region needs ad- 
ditional discussion. On one hand, it is required 
from mathematical point of view. Namely, the 
splitting of the velocity and density fields is ef- 
fective if all observed ions share the same veloc- 
ity distribution but respond differently to the gas 
density. If in addition to the varying density and 
velocity one would allow for the varying metallic- 
ity, the inverse problem becomes fully degenerate, 
i.e. it would have infinitely large number of solu- 
tions. On the other hand, the constant metallicity 
has some observational support: results obtained 
in numerous studies of the Galactic halo chem- 
ical composition reveal no systematic differences 
in the gas-phase abundances within galactocentric 
distances of 7 — 10 kpc in various directions (e.g., 
Savage & Sembach 1996). But, of course, we can- 
not exclude the case when the line of sight passes 
through many types of environments with differ- 
ent enrichment histories within a given absorber. 
If the metallicities within such an absorber differ 
only slightly ( < 0.5 dcx), the observed lines of 
different ions can be well fitted to the synthetic 
profiles calculated with some average value of the 
metal content. If, however, the differences in the 
metallicities are really large ( > 1 dex), the self- 
consistent fitting of all observed lines becomes im- 
possible. In this case we have to split the absorber 
into separate regions having different metal abun- 
dances (see § 4.2.2 and 4.3.1 for examples). 

It is well known that the measured metallici- 
ties depend in a crucial way on the adopted back- 
ground ionizing spectrum. We started in all cases 
with the Haardt-Madau (HM) background ioniz- 
ing spectra (Haardt & Madau 1996) computing 
the fractional ionizations and the kinetic temper- 
atures with the photoionization code CLOUDY 
(Ferland 1997). If the fitting with the HM spec- 
trum was impossible, we used another spectra, e.g. 
the Mathews and Ferland (MF) spectrum (Math- 
ews & Ferland 1987). 

The MCI procedure itself is implemented in the 
following way. Within the absorbing region the ra- 
dial velocity v(x) and the total hydrogen density 
nn(x) along the line of sight are considered as two 
random fields which are represented by their sam- 



pled values at equally spaced intervals Ax, i.e. by 
the vectors {v\, . . . , Vk} and {m, . . . , n^} with k 
large enough (<~ 150 — 200) to describe the nar- 
rowest components of the complex spectral lines. 
The radial velocity is assumed to be normally dis- 
tributed with the dispersion cr v , whereas the gas 
density is distributed log-normally with the mean 
n and the second central dimensionless moment 
cr y (y = nn/n ). Both stochastic fields are cal- 
culated using the Markovian processes (see LAK 
for mathematical basics). The model parameters 
estimated in the least-squares minimization of the 
objective function (see eqs.[29] and [30] in LAK) 
include cr v and <r y along with the total hydrogen 
column density Ah, the mean ionization param- 
eter Uq, and the metal abundances Z a for a ele- 
ments observed in a given absorption-line system. 

The computations are carried out in two 
steps: firstly a point in the parameter space 
(Ah, Uo, (t v , (T y , Z a ) is chosen at random and then 
an optimal configuration of {vi} and {iii} for this 
parameter set is searched for. These steps are re- 
peated till a minimum of the objective function 
(x 2 ~ 1 P er degree of freedom) is achieved. To 
optimize the configurations of {vi} and {n^}, the 
simulated annealing algorithm with Tsallis accep- 
tance rule (Xiang et al., 1997) and an adaptive 
annealing temperature choice is used (details are 
given in Paper I). 

The following important fact should be taken 
into account when one interprets the results ob- 
tained with the MCI technique. The mean ioniza- 
tion parameter Uq is related to the parameters of 
the gas cloud as (see eq.[28] in LAK) 

U =**.(l + o*). (1) 
n 

Here n p h is the number density of photons with 
energies above 1 Ry which is determined by the 
intensity of the adopted background ionizing spec- 
trum. 

This equation shows that if the density field 
is fluctuating (cr y > 0), then with the same mean 
density no and the same background ionizing spec- 
trum we obtain a higher value of Uo without any 
additional sources of the ionization. Intermittent 
regions of low and high ionization caused by the 
density fluctuations will occur in this case along 
the sightlinc. On the other hand, for a given Uo 
the mean gas density n is also higher in the fluc- 
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tuating media as compared to the completely ho- 
mogeneous gas clouds (cr y = 0). Since the linear 
size of the absorber is L = Ah /no, the sizes esti- 
mated with the assumption of a constant density 
(as, e.g., in the VPF) may occur too large. 

Another important question is whether the 
MCI solution is unique and accurate. In general, 
the inverse problems are highly non-linear and 
ill-posed which implies multiple solutions and/or 
very broad uncertainty ranges for the recovered 
parameters. To produce a physically reasonable 
solution, one has to account for all available in- 
formation related to the case under study. For 
instance, the more lines of different ionic tran- 
sitions are included into the analysis, the more 
accurate result can be obtained since both low- 
and high-density regions are probed simultane- 
ously. One may also compare the relative metal 
abundances predicted by nucleosynthctic theories 
with those provided by the MCI. An odd pattern 
may indicate a misleading solution. The recov- 
ered linear sizes must also be in agreement with 
the characteristic sizes of the absorbers stemming 
from observations of gravitationally lensed quasars 
and quasar pairs which show L < 100 kpc. 

One of the main problem in the analysis of 
the QSO high rcdshift spectra is the line blend- 
ing which hampers significantly the inversion of 
the observed spectra. As compared to the VPF 
method, the MCI is much more robust dealing 
with blended lines due to the assumption that all 
ions trace the same underlying gas density and ve- 
locity distributions. This means that we are able 
to reconstruct both distributions using unblended 
parts of different lines. It is obvious that the ac- 
curacy of the recovered parameters improves with 
increasing number of lines and the variety of ions 
involved in the analysis. A priori we do not know 
which parts of the lines are blended and which are 
not. To clarify this, several test runs with dif- 
ferent arrangements of lines are carried out till a 
self-consistent fit for the majority of lines observed 
in the spectrum is found. 

4. Results on individual metal systems 

All results given below in Tables 1-3 were ob- 
tained using the MCI procedure as described in 
Paper I. Giving the shape and the intensity at 1 
Ry of the background ionizing radiation, the errors 



of the fitting parameters Uq, Ah, cr v , a y and Z a are 
about 15-20% , the errors of the estimated column 
densities are less than 10%, whereas the derived 
parameters no and L are estimated with about 
50% accuracy. These errors, however, should be 
considered as internal in the sense that they reflect 
merely the configuration of the parameter space in 
the vicinity of a minimum of the objective func- 
tion. To what extent the recovered parameters 
may correspond to their real values is discussed 
separately for each individual absorbing system. 
We note in passing that the density n and by 
this the linear size L scales with the intensity of 
the radiation field (see eq.[l]). 

The analyzed metal systems arc described 
within three categories : (1) Lyman limit systems 
with A(H I) > 5 x 10 16 cm" 2 , (2) Lya absorbers 
with N(U I) < 5 x 10 16 cm~ 2 , and (3) Lya sys- 
tems with a probable metallicity gradient. Their 
physical properties are compared with different 
types of absorbers listed in Appendix. 

4.1. Lyman limit systems 

4.1.1. Q0347-3819, z abs = 2.21 

This system consists of a broad saturated Lya 
hydrogen line spread over 500 km s _1 and of 
metal lines of low and high ionized species : 
O IA1302, C IIA1334, Mg IIAA2796, 2803, 
Al IIA1670, Si IIAA1190, 1193, Fe IIAA1608, 2344, 
2382, 2586, 2600, Al IIIA1854, Si IIIA1206 and 
C IV and Si IV doublets as well. The physical pa- 
rameters obtained with the MCI are presented in 
Table 1 whereas the corresponding observed and 
synthetic spectra are shown in Fig. 2. The recov- 
ered density and velocity distributions along the 
line of sight are plotted in Fig. 3. The intermit- 
tent high- and low-density regions giving rise to, 
respectively, low- and high-ionization species (a 
multiphase medium) are clearly seen. This means 
that the lines of different ionization stages arise 
in different areas despite of having the same ra- 
dial velocities [see, e.g., the regions with v(x) ~ 
km s _1 in Fig. 3]. 

According to these data, the system under 
study is a compact (380 pc) warm cloud (Tki n — 
9000 K) with a high metal content (0.6 solar value 
for O) and a high velocity dispersion (cr v ~ 80 
km s _1 ). 

The inferred H I column density of 4.6 x 10 17 
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cm~ 2 classifies this system as a typical LLS [since 
log N(R I) > 17 cm -2 ]. In principle, this value lies 
beyond the applicability limit of the MCI which is 
formally valid for t 912 < 1 [if N(H I) = 4.6 x 10 17 
cm~ 2 , then t 9 i 2 = 3]. Besides, having only one 
saturated Lya line we cannot in any case say for 
sure that the estimated value is the real hydrogen 
column density. But there are also reasons that 
make the obtained solution quite plausible. 

First, we consider an absorber as a clumpy re- 
gion which implies that the ionizing radiation may 
penetrate the cloud from different directions with- 
out being significantly altered (i.e. we assume that 
the density of the background ionizing radiation is 
not reduced much and its spectral distribution is 
not changed considerably in a gas cloud which is 
not a slab of a uniform density). Second, the ob- 
served mixture of the O I and Fe II lines and the 
C IV and Si IV lines makes it possible to fix the 
mean ionization parameter [7 quite strictly be- 
cause the fractional ionization curves for, e.g., O I 
and C IV are very different. Thus, the H I column 
density can hardly be less than the value estimated 
by the MCI since in that case the metallicity would 
be higher than solar. Besides, the same velocity 
interval covered by both the Lya and Si III A1206 
lines suggests that we observe probably the real 
profile of the H I Lya (a blend with the Ly/3 line 
from the system at z a b s = 2.8102 which falls in the 
range —40 km s _1 < Av < 140 km s _1 has a lit- 
tle influence on the Lya red wing - see Fig. 2). 
Higher values of A(H I) cannot be excluded, but 
additional calculations have shown that the max- 
imum of the H I column density is limited : both 
available wings of the Lya line do not allow an 
increase in N(H I) by more than a factor of three. 

This uncertainty does not change, fortunately, 
the main characteristic of this system - we do ob- 
serve at z a b s = 2.21398 a compact (L < 1 kpc) 
metal-rich (Z > 0.1 Z©) cloud. As seen from 
Fig. 2, the synthetic profiles represent well all un- 
blended spectral features (broad absorption fea- 
tures at the position of the C IV doublet are not 
consistent with each other and, hence, they cannot 
be attributed to the real C IV profiles; the same is 
valid for the Si IV doublet). The measured relative 
abundances are discussed later in § 5.3. Here we 
note that a set of the identified species and their 
pattern as well as the estimated mass (~ 1O 4 M ) 
resemble parameters measured in so-called high 



and intermediate velocity clouds (HVC, IVC) in 
the Local Group (Wakker 2001). Although the 
nature of the HVCs and IVCs is still poorly un- 
derstood, they are unlikely to be phenomena iso- 
lated to the Milky Way. The system at z a b s =2.21 
may be one of similar HVCs likely to be encoun- 
tered in an high-rcdshift galactic halo. A more 
precise determination of its properties would re- 
quire the higher Lyman series lines to be included 
in the analysis (these lines, however, can be ob- 
served with space telescopes only). 

4.1.2. Q0347-3819, z abs = 2.81 

Here we identified three neutral hydrogen lines 
and several metal absorption lines both of low and 
high ionic transitions (see Fig. 4, the C IV doublet 
falls, unfortunately, in a wavelength coverage gap 
of the Keck spectrum). The estimated physical 
parameters are listed in Table 1 whereas the cor- 
responding synthetic spectra are shown in Fig. 4 
by the smooth curves. It should be noted, how- 
ever, that the solution with the H I column den- 
sity of 5.5 x 10 16 cm -2 is not unique because all 
available hydrogen lines are saturated and partly 
blended in the red wings. We found also another 
solution with N(H I) = 10 17 cm~ 2 . The solution 
presented in Table 1 has been chosen because it 
delivered a more or less consistent set of all pa- 
rameters: the velocity dispersion a v ~ 60 km s^ 1 , 
the mean gas density no ~ 10~ 3 cm~ 3 , the linear 
size L ~ 30 kpc, and the metallicity Z ~ 1/10Z©. 
For comparison, the second solution with 7V(H I) 
= 10 17 cm" 2 gives the metallicity Z ~ 1/20Z©, 
the size L ~ 70 kpc and the same other parame- 
ters. According to the both sets of the estimated 
parameters and the fact that rather strong low- 
ionization lines of C II and Si II are observed, the 
z a b s = 2.81 system could be related to an inner 
galactic halo. 

4.1.3. APM BR J0307-4945, z abs = 4-21 

In the available spectral range we identified hy- 
drogen lines Lya, Ly/3 and Ly7 as well as metal 
lines C IIA1334, Si IIA1526, C IIIA977, N IIIA989, 
Al IIIA1854, Si IIIA1206 and the C IV and Si IV 
doublets. From the metal line profiles, two main 
subsystems can be clearly distinguished: the first 
at v ~ —200 km s _1 , and the second at v ~ 200 
km s _1 (see Fig. 5). The analysis of the z a b s = 
4.21 absorber has been carried out in two steps: 
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(i) the subsystem at v ~ 200 km s -1 was treated 
separately (the results are listed in column 4 in 
Table 1 and the corresponding synthetic spectra 
are shown by the solid lines in Fig. 5), (ii) both 
systems were fitted together (column 5 in Table 1, 
the dotted lines in Fig. 5). 

The recovered metallicities are high (about 
1/3Zq at v ~ 200 km s _1 , and slightly higher 
at v ~ —200 km s _1 ) and their pattern is nearly 
solar. At low redshifts similar characteristics are 
measured, e.g., in high-mctallicity blue compact 
galaxies. At high redshifts, some radio galax- 
ies are known which show two spatially resolved 
emitting regions, slightly sub-solar metallicities 
and sizes of about tens of kpc, e.g., MRC 2104- 
242 at z = 2.49 (Overzier et al. 2001) or TN 
J1338-1942 at z = 4.11 (Dc Brcuck ct al. 1999). 
It is also known that galaxies showing morpholog- 
ical evidence of a merger have excessive velocity 
widths of their spectral lines (e.g. Mallen-Ornelas 
ct al. 1999). 

Thus we conclude that the LLS under study 
probably arises when the line of sight intersects 
two clumps which may be merging. The observed 
range of the metal lines (from ~ —300 to 400 
km s _1 ) and their very complex structures are 
also in line with this picture. 

Although the recovered values are self-consistent, 
we cannot guarantee their uniqueness because of 
the following reasons. The hydrogen lines are 
blended and saturated and, hence, the real N(R I) 
value may be higher leading to lower metal abun- 
dances. Besides, in our calculations the HM back- 
ground ionizing spectrum was adopted. However, 
it is probable that in close pair of galaxies where 
a significant star-forming activity is triggered by 
the merging, the HM spectrum can be affected by 
the local sources. In particular, the discrepancies 
between the observed and theoretical intensities 
seen in the C IV and Si IV components at v = 100 
km s _1 as well as the estimated overabundant 
ratio of [C/Si] ~ 0.1 (usually [Si/C] lies between 
and 0.4, see Appendix) may be caused by inade- 
quate choice of the background ionizing spectrum. 
Nevertheless, the observed metal profiles are quite 
consistent with the HM spectrum and therefore we 
expect that the influence of the local sources may 
not be very strong and, hence, the interpretation 
of this absorption system will not be significantly 
altered. 



4.2. Absorbers with A(H I) < 5 x 10 16 cm 2 

4.2.1. Q0347-3819, z ahs = 2.53 and 2.65 

These two metal systems with z a b s = 2.5370 and 
2.65044 show broad and saturated hydrogen Lya 
lines and lines of silicon and carbon in different 
ionization stages. Computational results are pre- 
sented in Table 2, the observed and synthetic pro- 
files are shown in Figs. 6 and 7. 

The solutions listed in Table 2 are non unique 
because of the sole hydrogen line and blended low- 
ionization lines of Si IIA1260 and C IIA1334. If 
in reality the Si IIA1260 absorption is absent at 
z a bs = 2.5370, then the solution with a higher Uq 
and, hence, with a large linear size can also be ob- 
tained. Therefore we consider the estimated lin- 
ear size of 13 kpc as a lower limit. The system at 
z a bs = 2.650 may have low-ionization C IIA1334 
absorption resulting in a lower mean ionization pa- 
rameter and a smaller linear size as compared with 
those listed in Table 2. A rather high overabun- 
dance ratio [Si/C] ~ 0.5 estimated for this system 
also allows us to speculate that the real ionization 
parameter may be lower. But these considerations 
do not change the classification of both absorbers: 
most probably they are hosted by halos of some 
distant galaxies. 

4.2.2. Q0347-3819, z abs = 2.962 and 2.966 

These two systems are separated by only 300 
km s _1 but demonstrate very different physical 
characteristics (see Table 2 and Figs. 8 and 9). 
Multiple hydrogen lines available in their spectra 
make it possible to estimate quite accurately the 
hydrogen column densities and all other physical 
parameters. 

The z a b s = 2.96171 system shows a broad Lya 
line extending over 400 km s _1 and weak absorp- 
tion lines of highly ionized silicon and carbon [the 
C IIIA977 line is contaminated by the Ly7 line 
from the z a b s = 2.97915 system (see Fig. 11); the 
C IIA1334 line is in a wavelength coverage gap of 
the Keck spectrum; the Si II AA 1260, 1193 lines 
are strongly blended]. The derived physical pa- 
rameters are typical for a galactic halo absorber: 
a low density (no ~ 3 x 10~ 4 cm~ 3 ), low metallic- 
ity (Z ~ 0.01 Z Q ), hot (T kin ~ 40000 K) cloud of 
L ~ 20 kpc size. The large overabundance of sili- 
con as compared with carbon ([Si/C] < 0.8) can 
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be explained by the uncertainty in the estimated 
carbon abundance: only one weak C IVA1548 line 
is available for the analysis. 

The adjacent system at z a b s = 2.96591 is on 
the contrary a very compact (L ~ 120 pc), warm 
(Tkin ~ 10000 K) cloud, 5 times more denser and 
30 times more metal abundant. This absorber re- 
veals also a weak N IIIA989 line [contaminated 
in the right wing by the H2 L11-0P(2) line from 
the z a bs = 3.025 DLA system - see Fig. 5 in Lev- 
shakov et al. 2002b]. We do not detect high 
amplitude fluctuations in the density and veloc- 
ity fields in this compact system (see Fig. fO), and 
as a result the observed metal line profiles are al- 
most symmetric. The low density region with the 
space coordinate < x < 0.2 does not contribute 
much to the line profiles, although a weak absorp- 
tion arising in this gas can be seen in Fig. 9 in 
the C IIIA977 and C IVA1548 lines at 11 ~ -30 
km s _1 . 

Since it is hardly possible that a cloud with a 
120 pc size could exist in space on its own, the 
two systems are probably physically related. For 
instance, this small and metal enriched cloud may 
be a condensation of a supernova-heated gas in a 
galactic halo seen in the z a b s = 2.96171 absorption 
lines. This process known as a galactic fountain 
(Brcgman 1980) is believed to give origin to the 
high mctallicity HVCs observed in the halo of the 
Milky Way. The velocity excess of the HVCs is 
usually greater than 90 km s _1 which is consis- 
tent with the shift of ~ 120 km s _1 between the 
redward absorption in the Lya profile at z a b s = 
2.96171 and the centre of the Lya line at z a b s = 
2.96591. 

4.2.3. Q0347-3819, z abs = 2.98 

Although this system belongs to the most com- 
monly observed type in the Lya forest (apart from 
hydrogen lines only weak C IV and Si IV doublets 
and no apparent absorption in other ionic species 
are registered), its hydrogen lines are unusually 
broad - extended over 600 km s _1 . 

Available Lya, Ly/3 and Ly7 lines (see Fig. 11) 
allow us to estimate the hydrogen column density 
with a sufficiently high accuracy. According to 
the recovered physical parameters (see Table 2), 
in this case we are dealing with a large (L ~ 50 
kpc) cloud of a rarefied (n — 10~ 4 cm~ 3 ), metal- 



poor (Z < 0.01 Z e ) and hot (T kin ~ 40000 K) gas. 
Blending of the C IIIA977 and Si III A1206 lines 
and the weakness of the Si IV A 1393 line do not 
allow us to estimate accurately the mean ioniza- 
tion parameter. The Uq value presented in Table 
2 should be considered as a lower limit. If in re- 
ality Uo is higher, then the absorber may have a 
lower mean gas density and a larger linear size. 

The most probable host for this system might 
be an external region of a giant galactic halo or a 
large scale structure object. 

4.24. Q0347-3819, z abs = 3.14 

The unsaturated Ly7 line gives accurate esti- 
mations of the total neutral hydrogen column den- 
sity (Fig. 12, Table 2). Clean continuum windows 
seen at the expected positions of metal lines make 
it possible to estimate the upper limits on metal 
abundances and to calculate the total hydrogen 
column density TVjj. The result obtained shows a 
rather low mctallicity cloud with [C/H] < —2.2 
and the linear size L > 13 kpc. One may expect 
to observe similar systems in the outer parts of 
galactic halos. 

4.2.5. APM BR J0307-4945, z abs = 4. 81 

This is the most distant absorber in our set 
where a low metal abundance can be directly mea- 
sured. Its hydrogen Lya line is clearly seen at 
Av = 3000 km s _1 in the wide emission blend 
Lya+N V+Si II shown in Fig. lc. Absorption 
lines in Fig. 13 give the redshift z a b s = 4.8101 and, 
thus, this system has z a b s > z om . The same or- 
der of magnitude velocity difference (Av ~ 3000 
km s _1 ) between the H 2 -bearing cloud at z a b s = 
2.811 (Levshakov & Varshalovich 1985) and the 
quasar redshift z cm = 2.770 (Foltz et al. 1988) 
has been observed toward PKS 0528-250. This H 2 
cloud seems to be at a distance larger than 10 kpc 
from the quasar as shown by Srianand & Petit- 
jean (1998). However, in our case we are not able 
to estimate the proximity of the quasar. We can 
only assume that the photoionization of the z a bs = 
4.8101 system could be affected by the quasar ra- 
diation. 

This assumption is supported by the follow- 
ing facts. It turned out to be impossible to 
fit all available lines with the HM ionizing spec- 
trum: the relative intensities of the C IIIA977 
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and C IVA1548,1550 required very high Uq val- 
ues (Uo > 0.1) which contradicted with the shal- 
low extended wings of the Lya line. On the con- 
trary, the MF ionizing spectrum, corresponding to 
the AGN, allowed to fit all lines and delivered a 
self-consistent set of physical parameters (Table 2, 
Figs. 13 and 14). The intensity at 1 Ry was set 
to J912 = 1CP 22 erg cm~ 2 s _1 Hz -1 sr _1 which 
corresponds to the intensity of the HM spectrum 
at z = 4.9. 

According to the recovered values, the absorber 
at z a b s = 4.81 is a metal-poor cloud with the linear 
size of about 25 kpc and the mean density tiq ~ 2 x 
1CP 4 cm~ 3 . Fig. 14 shows that the shallow wings 
of the Lya line are produced by the streaming out 
low density gas whereas the central region of the 
cloud remains very quiet. 

Since we cleary see only the C IV doublet 8 and 
have upper limits for the intensity of the C III A977 
and Si IV A1393 lines, the mean ionization param- 
eter Uo listed in Table 2 should be considered as a 
lower limit implying that lower metal abundances 
and lower gas densities may also be possible. Tak- 
ing into account that the density no and by this 
the linear size L scales with the intensity of the 
radiation field (see eq.[l]), the value of L becomes 
quite uncertain. It could be larger or smaller than 
the estimated size of 25 kpc (this uncertainty is 
marked by '?' in Table 2). 

An absorption system with similar spectral 
characteristics (wide and shallow wings of the 
Lya line, a weak C IV doublet) was previously 
observed by Reimers et al. (2001) at z a b s = 1.674 
toward HE 0515-4414 (z om = 1.73). The system 
at z a b s = 1-674 shows in addition a strong O VI 
doublet which, unfortunately, cannot be identified 
in the z a bs = 4.81 cloud because of blending in the 
dense Lya forest. 

Thus we may conclude that such systems are 
probably formed in the gas clouds affected by the 
QSO radiation. 

The A1550 A component is slightly blended with the night 
sky lines, whereas the A1548 A line is clear as found from 
the comparison of the two spectra of the APM BR J0307- 
4945 taken in a two month interval. 



4.3. Lya systems with metallicity gradient 

4.3.1. Q0347-3819, z abs = 2.848 and 2.899 

These systems present several seemingly un- 
blended hydrogen lines and pronounced lines 
of highly ionized silicon and carbon (Figs. 15, 
16; Table 3). Additionally, the N VA1238 and 
O VI A 1031 lines were identified at z abs = 2.848. 

In spite of the multiple H I lines, it turned out 
that the MCI failed to fit adequately all available 
ionic profiles in the apparent velocity ranges when 
homogeneous metallicities over the entire absorb- 
ing regions were assumed. The most sensitive re- 
strictions in these calculations are set by the con- 
tinuum windows seen in the profiles of the strong 
C IIIA977 and Si IIIA1206 lines. We consider 
the observed lines from these systems as arising 
in clouds with metallicity gradients for the follow- 
ing reasons. 

The hydrogen Lya and Ly/3 lines seen at z a b s = 
2.848 in the Av range between 100 km s _1 and 
400 km s _1 do not show any metal absorption 
in this velocity interval. Absorption features seen 
in the range 140 km s _1 < Av < 400 km s _1 
in Fig. 15 (panel C III) and in the range 180 
km s _1 < Av < 400 km s _1 in panel Si III cannot 
be attributed to the corresponding hydrogen sub- 
system, otherwise we should observe pronounced 
C IV lines in the same velocity range. If the red- 
ward hydrogen absorption is physically connected 
with the blueward one, then we may conclude that 
the z a bs = 2.848 absorber has a metallicity gradi- 
ent. 

A very similar picture is seen in the z a bs = 2.899 
system for the blueward portions of the hydrogen 
lines in the range —220 km s _1 < Av < — 80 
km s- 1 (Fig. 16). The C III A977 line (which is 
the most sensitive absorption line in a wide range 
of the ionization parameter from log Uo > — 4 to 
log [/0 — 0) shows no absorption in this velocity 
range. 

In principle, the z a b s = 2.899 system may be 
caused by a blending effect when the line of sight 
intersects two separate clouds (one of them sup- 
posingly metal-free) which have approximately the 
same radial velocities. As for the z abs = 2.848 sys- 
tem, a fortuitous blending can be ruled out be- 
cause we observe two absorption features of ap- 
proximately equal depth at v ~ —80 km s _1 in 
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the Lya and C IV profiles - the configuration no- 
ticed previously at z a bs = 1.385 toward HE 0515- 
4414 by Reimers et al. (2001), i.e. this unusual 
pattern is not unique and seems to indicate some 
special kind of absorption systems. Preliminarily 
we may conclude that the system at z a b s = 2.848 
as well as probably that at z a b s = 2.899 arise when 
the line of sight intersects a distant halo with very 
low (if any) metal content and encounters a metal- 
rich HVC. Unfortunately, the accurate quantita- 
tive analysis of both systems cannot be carried 
out because the velocity excesses of the assumed 
HVCs are not large enough to separate the inputs 
to the hydrogen lines from the halos and from the 
clouds like it was possible in the case of the z a bs = 
2.962 and z a b s = 2.965 systems. Therefore the re- 
sults presented in Table 3 are tentative, but they 
are physically reasonable and self-consistent. 

The corresponding synthetic spectra are shown 
in Figs. 15 and 16 by the solid lines. It should 
be noted that for both systems, only one wing of 
the available H I lines (marked with the horizontal 
bold lines) was included in the analysis, whereas 
the synthetic profiles for the entire H I lines of 
the encountered HVCs were computed using the 
velocity and density distributions estimated from 
the metal profiles. 

According to the data from Table 3, the sug- 
gested HVCs belong to different types. The sys- 
tem at z a bs = 2.899 has the metal abundance pat- 
tern, the set of the ions observed (notice the 
pronouced Si IV doublet) and the size similar to 
those estimated for the supposed HVC at z a b s = 
2.965 and consistent with the parameters of the 
HVCs observed in the Milky Way (see Appendix) . 
The system at z a b s = 2.848 is more higher ionized 
(O VI and N V and no Si IV) and it has the size 
of several kiloparsecs. Highly ionized HVCs with 
similar parameters were observed by Sembach et 
al. (1999) near the Milky Way. Their origin is 
still uncertain. Probably the z a b s = 2.848 absorber 
may belong to the intercluster gas clouds in a dis- 
tant group of galaxies as it was suggested by Sem- 
bach et al. for the local highly ionized HVCs. 

5. Discussion 

5.1. The origin of metal systems 

Metal systems with N(H I) < 5x 10 16 cm" 2 are 
usually believed to originate in galactic halos at 



different galactocentric distances. At low redshifts 
(z < 1) the galaxies associated with certain metal- 
lic absorptions (e.g. C IV) can be in most cases 
identified directly (e.g. Chen et al. 2001a). Our 
results on absorption systems with z > 2 also sup- 
port this assumption: absorbers with z a b s = 1.87 
(Paper I), 2.54, 2.65, 2.962, 2.98 (present paper) 
are produced by metal-enriched [Z < 0.1 Zq), hot 
(Tkin > 20000 K), rarefied (n ^ 10~ 4 - 10~ 3 
cm~ 3 ) gas clouds which have typical linear sizes 
of L > 10 kpc. These parameters are consistent 
with contemporary models of galactic halos (e.g. 
Viegas, Friaca & Gruenwald 1999). 

The nature of Lyman limit absorbers is less 
understood. Mo & Miralda-Escude (1996) as- 
sociate them with cold photoionized clouds ran- 
domly moving in hot spherical halos. The clouds 
are supposed to form from the initial density inho- 
mogeneities in the accreting intergalactic gas dur- 
ing its cooling. Both the cloud and the halo ob- 
viously reveal the equal metallicity since they are 
formed from the same gas. In our study two LLSs 
with z a b s = 1-92 (Paper I) and 2.81 (present paper) 
can be related to the absorbers of this type. 

However, this scenario obviously fails to explain 
metal abundant (Z > 0.1 Zq) systems since it is 
hard to understand how the whole halo can be 
metal-enriched to such a high level. It was shown 
by hydrodynamic simulations (e.g. Katz et al. 
1996; Gardner et al. 2001) that LLSs can also 
arise on lines of sight that pass through small pro- 
togalaxies. We found two systems with z a bs = 1-94 
(Paper I) and 4.21 (present paper) that can be ex- 
plained within this framework. These metal-rich 
(Z ~ 1/3 Zq) absorbers with the sizes of several 
kpc are probably hosted by objects that may be 
akin to the local compact blue galaxies. 

Some absorbers in our present study (z a b s = 
2.21, 2.965, and, possibly, 2.89) reveal small lin- 
ear sizes (L < 1 kpc) together with very high 
metal content (Z ~ 1/2 Zq). These three sys- 
tems may be explained in the framework of the 
process known as a galactic fountain : metal- 
enriched (supernova-heated) gas arises from the 
inner region of a galaxy and condenses into the 
clouds within the hot galactic halo. After for- 
mation, clouds fall back toward the galaxy cen- 
tre because of their higher density. It is supposed 
that high-metallicity HVCs observed in the Milky 
Way halo are formed by this mechanism (Breg- 
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man 1980). The HVCs are common objects in our 
Galaxy and are detected in every longitude and 
latitude region. If galactic fountain works also in 
distant galaxies, it would be quite probable to en- 
counter such a cloud on the line of sight which 
intersects the galactic halo, as also discussed by 
Charlton et al. (2001). Another type of HVCs - 
hot, highly ionized clouds with sizes of several kilo- 
parsecs - is represented by the absorption system 
at z a b s = 2.848. The origin of this type of HVCs is 
uncertain, but they may be produced by the inter- 
galactic metal-enriched gas falling onto metal-poor 
galactic halos. 

Measured abundances of C and Si are depicted 
versus logarithmic sizes of the studied systems in 
Fig. 17. Systematically higher metal abundances 
are seen in compact systems with linear sizes L < 
4 kpc. This result seems to indicate that the more 
effective metal enrichment occurs within relatively 
compact regions. 

Our results show that Lyman limit systems are 
a heterogeneous population which is formed in at 
least three different environments. This should be 
taken into account when statistics of LLSs is used 
to verify different models in hydrodynamic cosmo- 
logical simulations. 

5.2. (7 V — A H L relation 

If QSO metal systems are formed in gas clouds 
gravitationally bound with intervening galaxies, 
the internal kinematics of the QSO absorbers 
should be closely related to the total masses of 
the host galaxies. In case of galactic population, 
different types of galaxies show different scaling 
relations between the linear size and the velocity 
width of emission lines (e.g., Mallen-Ornelas et 
al. 1999). Possible correlation between the ab- 
sorber linear size L and its line-of-sight velocity 
dispersion cr v was also mentioned in Paper I. 

The correlation between o~ v and L stems from 
the virial theorem which states : 

<r*~ ^-~n L 2 =NnL. (2) 

Assuming that the gas systems are in quasi- 
cquilibrium, one can expect cr v ~ V^h L. 

In Fig. 18 we examine our systems by compar- 
ing their kinematics (ct v ) with measured sizes (L) 
and total gas column densities (-/V H )- Shown are 
the data for all QSO absorbers studied in Paper I 



and in the present paper except for the systems 
at z abs = 2.848 and 2.899 (Table 3) which show 
inhomogeneous metallicities. It is seen that in the 
log(er v ) versus log(A^H L) diagram, most systems 
with linear sizes L > 1 kpc lie along the line with 
the slope k = 0.30 ± 0.03 (1 a c.L). 

Taking into account that we know neither the 
impact parameters nor the halo density distribu- 
tions, this result can be considered as a quite good 
fit to the relation (1). Hence we may conclude that 
most absorbers with L > 1 kpc are gravitationally 
bound with systems that appear to be in virial 
equilibrium at the cosmic time when the corre- 
sponding Lya absorbers were formed. The possi- 
ble consequence of this conclusion is that since the 
most metal rich absorbers identified in the QSO 
spectra arise in the galactic systems the question 
whether the intergalactic matter is metal enriched 
or pristine remains still open. 

5.3. [a-element/iron-peak] ratio 

The metal abundances measured in the z a b s = 
2.21 LLS (Table 1) can be used to estimate the 
a-element to the iron-peak group ratio which is a 
good indicator of the chemical evolutionary status 
of high rcdshift gas clouds. During the chemical 
evolution, heavy elements produced in stars show 
different nucleosynthetic histories so that their rel- 
ative abundances vary with cosmic time. 

Oxygen and other a-chain elements are mainly 
produced by Type II SNe, while iron is also a 
product of Type la SNe which have longer evo- 
lution scales. In the early stages of the chemical 
evolution of galaxies (At < 2 x 10 7 yr) the in- 
terstellar gas is likely enriched by Type II SNe 
products, while at At > 10 8 yr, the [a/Fe] ra- 
tio should decline. Observations reveal both low 
[e.g. ~ 0.1 — 0.2 in the z a b s = 3.390 dust-free 
DLA (Q0000-2621; Molaro et al. 2001) and in 
the z abs = 3.386 DLA (Q0201+1120; Ellison et 
al. 2000)], and high [e.g. ~ 0.7 in the DLA 
I Zw 18 (Levshakov, Kegel & Agafonova 2001) and 
0.68 ± 0.08 in the z abs = 3.025 DLA (Q0347-3819; 
Levshakov et al. 2002b)] ratios of [a-element /iron- 
peak] . 

Oxygen with its weak affinity with dust grains 
is a good tracer of the a-element abundances. 
Nevertheless, the intrinsic [a/Fe] ratio may be af- 
fected by depletion of iron since being a refrac- 
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tory element iron may be partly locked into dust 
grains. The dust content in the z a bs = 2.21 LLS 
may not, however, be too high. The relative abun- 
dances of the a-elements O, Mg and Si are [Si/O] 
= -0.14 ±0.11 and [Mg/O] = -0.27 ±0.11. In 
Galactic stars the a-elements show the same be- 
haviour relative to iron-peak elements (ovcrsolar 
at [Fe/H] < — 1; see, e.g., Gaswami & Prant- 
zos 2000) . We thus expect to find solar a-element 
ratios in dust-free absorbing regions, as observed, 
e.g., in the mentioned above z a bs — 3.390 DLA to- 
ward Q0000-2620. A negative value of [Mg/O] 
found in this LLS may indicate the presence of 
some amount of dust with a depletion factor of 
about 0.2 dex for the magnesium abundance. If, 
however, only the gas-phase abundances of O and 
Fe are taken, the upper bound on the [O/Fe] ra- 
tio is 0.65 ± 0.11, which is comparable with that 
found, for instance, in the z a b s = 3.025 DLA to- 
ward Q0347-3819 where the dust-to-gas ratio is 
~ 1 /30 of the mean Galactic interstellar medium 
value (Levshakov et al. 2002b). The enrichment of 
the a-elements in the z a b s =2.21 LLS is also sup- 
ported by the relative abundances of Si, Mg to Fe: 
[Si/Fe] = 0.51 ±0.11 and [Mg/Fe] = 0.38 ±0.11. 
Thus, the absorbing cloud at z a b s = 2.21 appears 
to be a chemically young object. 

6. Summary 

We have deduced the physical properties of ten 
absorption- line systems in the range Az = 2.21 — 
2.966 toward Q0347-3819 and of two systems at 
z abs = 4.21 and 4.81 toward APM BR J0307-4945. 
The main conclusions are as follows : 

1. The analyzed Lyman limit systems belong to 
a heterogeneous population which is formed 
by at least three groups of absorbers : (i) ex- 
tended metal-poor gas halos of distant galax- 
ies; (ii) gas in dwarf galaxies; {Hi) metal- 
enriched gas arising from the inner galac- 
tic regions and condensing into the clouds 
within the hot galactic halo (galactic foun- 
tain). While the interpretation of a single 
system is sometimes subject to large uncer- 
tainties as discussed in chapter 4, the exis- 
tence of a wide spread of properties in the 
different systems is certainly proved. 

2. A correlation between the line-of-sight ve- 
locity dispersion ct v and the linear size L of 



the absorbing systems noted in Paper I is 
confirmed. New results show that large size 
QSO absorbers (L > 1 kpc) obey a scaling 
law a v ~ (iVn L) 3 over two decades in ve- 
locity dispersion and in the product (Ah L). 
This means that the majority of the metal 
absorbers are probably bound to the galac- 
tic systems and, hence, the question whether 
the IGM is enriched or pristine requires fur- 
ther investigations. 

3. Systematically higher metal abundances are 
found in compact systems: in our sample 
there are no small size systems (L < 1 kpc) 
with metallicity lower than 0.1 Z@. 

4. The gas-phase metal abundances from the 
z a bs = 2.21 LLS reveal a pronounced fa- 
element /iron-peak] enhancement with [O/Fe] 
= 0.65 ± 0.11 at the 6er confidence level, the 
first time when this abundance pattern is 
unambiguously found in a LLS. The mea- 
sured [O/Fe] ratio implies that the chemical 
history of this LLS is < 10 8 yr. 

5. Absorption system at z a b s = 2.21 and 2.965 
and possible the systems at z a b s = 2.848 
and 2.899 toward Q0347-3819 show char- 
acteristics very similar to that observed for 
different types of HVCs in the Milky Way 
and may be interpreted as being the high- 
redshift counterparts of these Galactic ob- 
jects. 
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edges the hospitality of the National Astronomical 
Observatory of Japan (Mitaka) where this work 
was performed. The work of S.A.L. and I. LA. is 
supported in part by the RFBR grant No. 00-02- 
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A. Possible counterparts of QSO absorbers 



In our interpretation of the nature of the different metal systems with iV(H I) ~ 10 14 — 10 17 cm -2 we have 
referred to a set of absorbers whose physical parameters are summarized in Table 4. Comparison between 
the reference absorbers and the QSO systems is based on following parameters : 

1. Linear sizes. The systems can be divided into small-size (L < 1 kpc), intermediate-size (1 kpc < L < 10 
kpc), and large-size (10 kpc < L < 150 kpc) absorbing regions connected with galactic gaseous 
envelopes, as well as into very large-size (L > 150 kpc) absorbers (filaments) showing correlation with 
the large-scale distribution of galaxies (e.g. Penton et al. 2002). Filaments are probably intergalactic 
material not recycled by galaxies. 

2. Metallicities. Values of 0.3 < Z/Z & < 1, 0.03 < Z/Z Q < 0.3, and Z/Z Q < 0.03 classify the systems as 
metal rich, metal enriched, and metal poor, respectively. 

3. Metallicity patterns. Relative element abundances allow us to distinguish between chemically young 
and old systems (low and high relative [a-element/iron-peak] ratio, respectively). 
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Fig. 1 — Emission line profiles in the VLT/UVES 
spectrum of APM BR J0307-4945. The zero ve- 
locity corresponds to the rcdshift z om = 4.7525 
and indicates the adopted line centre of the O I 
A1304 line. A smooth fit to the C IV A1549 pro- 
file superimposed on the original line is shown in 
panel (b) . Comparison of the synthetic C IV pro- 
file scaled to match other emission lines is illus- 
trated in panels (a), (c), (d), and (e). A combined 
profile constructed from the same synthetic C IV 
profile properly scaled for each emission line and 
shifted for the N V A1240 and Si II A1264 lines 
to their centers (assuming z em = 4.7525) is super- 
posed on the whole extended emission line in panel 
(*)■ 



Fig. 2.— Hydrogen and metal absorption lines as- 
sociated with the z a b s = 2.21 LLS toward Q0347- 
3819 (normalized intensities are shown by dots 
with la error bars). The zero radial velocity 
is fixed at z = 2.21398. Smooth lines are the 
synthetic spectra convolved with the correspond- 
ing point-spread functions and computed with the 
physical parameters from Table 1. Bold horizon- 
tal lines mark pixels included in the optimization 
procedure. Spectra obtained with the VLT/UVES 
and the Keck/HIRES spectrographs are marked 
by VLT and Keck, respectively The normalized 
Xmin = 1-39 (the number of degrees of freedom 
v = 1255). 



Fig. 3. — Computed velocity (upper panel) and 
gas density (lower panel) distributions along the 
line of sight within the z a bs = 

2.21 LLS toward 
Q0347-3819. Shown are patterns rearranged ac- 
cording to the principle of minimum entropy pro- 
duction rate (see Paper I). 



Fig. 4. — Same as Fig. 2 but for the z a b s = 
2.81 LLS. The zero radial velocity is fixed at 
z = 2.8102. The corresponding physical param- 
eters are listed in Table 1. Xmin = 1-38, v = 763. 



Fig. 5. — Same as Fig. 2 but for the z a b s = 4.21 
LLS toward APM BR J0307-4945. The zero ra- 
dial velocity is fixed at z = 4.2144. The physical 
parameters listed in Table 1, column 4, correspond 
to the system centered at v ~ 200 km s _1 . The 
synthetic spectra for this system are shown by the 
solid lines (x^m = !- 83 > v = 1436). The synthetic 
spectra shown by the dotted lines are calculated 
with the parameters listed in Table 1, column 5, 
which are obtained from the analysis of the entire 
system from v = —300 km s _1 to 400 km s _1 
(xLn = 2-61, v = 2325). 

Fig. 6. — Same as Fig. 2 but for the z a b s = 
2.53 system. The zero radial velocity is fixed at 
z = 2.5370. The corresponding physical parame- 
ters are listed in Table 2. Xmin = L07 : v = 602 - 



Fig. 7. — Same as Fig. 2 but for the z a b s = 2.65 
system. The zero radial velocity is fixed at z = 
2.65044. The corresponding physical parameters 
are listed in Table 2. Xmin = L15 > v = 962 - 

Fig. 8. — Same as Fig. 2 but for the z a b s = 2.962 
system. The zero radial velocity is fixed at z = 
2.96171. The corresponding physical parameters 
are listed in Table 2. Xmin = !- 20 > v = 627 - 

Fig. 9. — Same as Fig. 2 but for the z abs = 2.966 
system. The zero radial velocity is fixed at z = 
2.96591. The corresponding physical parameters 
are listed in Table 2. Xmin = L10 > v = 149 - 

Fig. 10. — Computed velocity (upper panel) and 
gas density (lower panel) distributions along the 
line of sight within the z a b s = 2.966 absorber to- 
ward Q0347-3819. Shown are patterns rearranged 
according to the principle of minimum entropy 
production rate (see Paper I). 

Fig. 11. — Same as Fig. 2 but for the z abs = 2.98 
system. The zero radial velocity is fixed at z = 
2.97915. The corresponding physical parameters 
are listed in Table 2. Xmin = 1 - 20 > v = 858 - 
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Fig. 12. — Same as Fig. 2 but for the z abs — 3.14 
system. The zero radial velocity is fixed at z = 
3.13985. The corresponding physical parameters 
are listed in Table 2. xL in = v = 407 - 



Fig. 13. — Same as Fig. 2 but for the z abs = 4.81 
Lya system toward APM BR J0307-4945. The 
zero radial velocity is fixed at z = 4.8101. The cor- 
responding physical parameters are listed in Ta- 
ble 2. xLn = 0-85, " = 303. 



Fig. 14. — Computed velocity (upper panel) and 
gas density (lower panel) distributions along the 
line of sight within the z a b s = 4.81 absorber to- 
ward APM BR J0307-4945. Shown are patterns 
rearranged according to the principle of minimum 
entropy production rate (see Paper I). 



Fig. 15. — Same as Fig. 2 but for the z a b s = 2.84 
Lya system with inhomogeneous metallicity. The 
zero radial velocity is fixed at z = 2.84829. The 
corresponding physical parameters are listed in 
Table 3. xlin = 1-29, v = 284. 



Fig. 16. — Same as Fig. 2 but for the z a b s = 2.90 
Lya system with inhomogeneous metallicity. The 
zero radial velocity is fixed at z = 2.89922. The 
corresponding physical parameters are listed in 
Table 3. Xmin = L35 > v = 82 0- 



Fig. 18. — Plot of the line of sight velocity dis- 
persion log(<r v ) vs. log(7VH L) for the sample of 
metal absorption-line systems. The dashed line is 
the linear regression log(cr v ) oc k log(7Vu L) calcu- 
lated for the points shown as filled circles having 
both horizontal error bars. The slope k is equal to 
0.30±0.03 (1 cr c.l.). Shown are the data from Ta- 
bic 1 (Paper I) and form Tables 1 and 2 (present 
study). The data for the z a b s = 4.21 system to- 
ward APM BR J0307-4945 are taken from column 
4, Tabic 1. The point at log (N K L) = 41.08 cor- 
responds to the z a b s = 4.81 system toward APM 
BR J0307-4945 where the linear size is very un- 
certain. Open squares represent HVC-like systems 
(see text for details). 



Fig. 17. — Abundances of carbon (circles) and 
silicon (squares) plotted against the logarithmic 
linear size of the corresponding systems as mea- 
sured by the MCI procedure. Shown are the data 
from Table 1 (Paper I) and form Tables 1 and 2 
(present study). The data for the z a b s = 4.21 sys- 
tem toward APM BR J0307-4945 are taken from 
column 4, Table 1. The point at log (L) = 3.40 cor- 
responds to the z abs = 4.81 system toward APM 
BR J0307-4945 where the linear size is very un- 
certain. Small size QSO absorbers seem to have 
systematically higher metal content than systems 
with larger sizes. 
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Table 1 

Physical parameters of the Lyman limit systems toward Q0347-3819 and APM BR 
J0307-4945 derived by the MCI procedure 





Q0347-3819 


Q0347-3819 


APM BR J0307-4945 


Parameter" 


z ab3 = 2.21398 


z ab3 = 2.8102 


2 abs — 


4.2144 


(i) 


(2) 


(3) 


(4) 


(5) 


U 


2.1E-3 


4.0E-2 


7.5E-3 


8.1E-3 


Nn, cm~ 2 


2.2E19 


8.3E19 


4 2E19 


9 1E19 


<t v , km s~ 1 


80.8 


59.5 


80.0 


170.0 


°"y 


1.0 


1.15 


1.2 


1.25 


7 b 
z c 


2.0E-4 


3.1E-5 


1.3E-4 


2.0E-4 


Z N 




<4.0E-6 


<2.5E-5 


<2.5E-5 


Zo 


4.4E-4 








Ztjig 


1.5E-5 








Zai 






<2.5E-6 


<2.5E-6 


Zsi 


2.0E-5 


4.2E-6 


1.1E-5 


1.5E-5 


Zfo 


5.0E-6 








[Zc] c 


-0.28 


-1.10 


-0.47 


-0.27 


[Zn] 




< -1.3 


< -0.5 


< -0.5 


[Zo] 


-0.10 








[Z Ms ] 


-0.37 








[Zai] 






< -0.5 


< -0.5 


[Z Si ] 


-0.24 


-0.90 


-0.51 


-0.36 


[ZFe] 


-0.75 








JV(Hi), cm- 2 


4.6E17 


5.5E16 


1.4E17 


2.3E17 


N(Oi) 


1.5E14 








Af(Cn) 


1.5E15 


5.2E13 


4.8E14 


1.2E15 


JV(Mgn) 


1.1E14 








JV(Sin) 


2.3E14 


8.2E12 


4.8E13 


1.0E14 


JV(Feii) 


1.3E13 








AT (Cm) 




1.6E15 


4.3E15 


1.3E16 


AT(Niii) 




<2.0E14 


<6.3E14 


<1.6E15 


iV(Alm) 


<8.1E12 




<1.1E13 


<3.4E13 


AT(Sim) 


2.0E14 


1.2cl4 


2.6E14 


6.8E14 


N(C iv) 


7.7E13 




5.6E14 


2.0E15 


AT(Siiv) 


3.0E13 


6.7E13 


9.0E13 


2.9E14 


no, cm -3 


2.0E-2 


9.2E-4 


2.5E-3 


2.5E-3 


(T>, K 


9.1E3 


2.4E4 


1.2E4 


1.2E4 


T min 


8.6E3 


1.7E4 


1.1E4 


1.1E4 


T max 


9.8E3 


3.3E4 


1.4E4 


1.5E4 


L, kpc 


0.38 


30 


6.0 


13 



Note. — "The internal errors of Uq, Nh, cr v , tr y and Z a arc c± 15 — 20%, the 
column density errors arc < 10%, whereas no and L are known with c± 50% 
accuracy; b — X/H; c [^x] — log(Zx) — logfZx)© [solar abundances are 
taken from Holwcgcr (2001) except AI for which solar abundance from Grcvesse 
& Sauval (1998) is used ]; 
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Table 2 

Physical parameters of systems with iV(Hi) < 5 x 10 16 cm~ 2 toward Q0347-3819 (A) and 
APM BR J0307-4945 (B) derived by the MCI procedure 





A, Z abs — 


A, z abB = 


A, z aba = 


A, z aba = 


A, z abB = 


A, z abs = 


B, z abB = 


Parameter" 


2.5370 


2.65044 


2.96171 


2.96591 


2.97915 


3.13985 


4.8101 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Ua 


2.9E-2 


4.2E-2 


8.5E-2 


2.0E-2 


0.15 


0.11 


3.8E-2 


Nu, cm~ 2 


4.4E19 


1.9E19 


2.0E19 


5.1E17 


3.6E19 


1.2E19 


1.6E19 


ct v , km s _1 


68.5 


73.0 


88.0 


11.6 


110.8 


25.3 


33.8 


°y 


1.05 


0.7 


0.85 


0.92 


1.12 


1.0 


0.87 


7 b 
z c 


5.0E-6 


3.7E-5 


2.0E-6 


1.5E-4 


3.1E-6 


< 1.8E-6 


3.3E-6 


Z N 








1.8E-5 








Z S i 


1.2E-6 


1.2E-5 


1.1E-6 


3.2E-5 


<3.5E-6 


< 9.0E-7 


< 7.0E-7 


[ZcV 


-1.89 


-1.02 


> -2.3 


-0.42 


-2.07 


< -2.3 


< -2.0 


[Zn] 








-0.67 








[Zsi] 


-1.45 


-0.47 


-1.49 


-0.03 


< -1.0 


< -1.6 


< -1.7 


JV(Hl), cm- 2 


2.2E16 


6.2E15 


2.0E15 


3.6E14 


1.9E15 


8.7E14 


3.7E15 


JV(Cn) 








9.4E11 




<4.7E10 


<1.5E11 


Af(Sin) 


<1.2E12 






<1.8E11 




<1.6E10 




AT(Ciii) 






1.1E13 


4.3E13 


1.8E13 


<4.6E12 


< 4.8E12 


JV(Nin) 








5.2E12 








JV(Sira) 


1.6E13 


2.7E13 


1.1E12 


3.6E12 


<9.1E11 


<2.8E11 




iV(Civ) 


4.7E13 


2.5E14 


8.4E12 


1.9E13 


2.0E13 


<4.5E12 


9.4E12 


iV(Siiv) 


7.3E12 


3.7E13 


8.5E11 


5.2E12 


<1.4E12 


<3.0E11 


< 6.7E11 


no, cm~ 3 


<l.lE-3 


5.7E-4 


<3.2E-4 


1.5E-3 


<2.4E-4 


<3.1E-4 


2.3E-4? 


<T>, K 


2.8E4 


1.8E4 


4.0E4 


1.0E4 


4.4E4 


3.4E4 


3.2E4 


T min 


1.9E4 


1.4E4 


2.6E4 


7.0E3 


3.0E4 


2.7E4 


2.2E4 


^max 


4.3E4 


2.3E4 


5.0E4 


2.1E4 


7.0E4 


5.3E4 


4.3E4 


L, kpc 


>13 


11 


>20 


0.12 


>50 


>13 


25? 



NOTE. — °Thc internal errors of Uo, Nh, ct v , cry and Z a arc ~ 15 — 20% , the column density errors are < 10%, 
whereas no and L are known with ~ 50% accuracy; Zy^ — X/H; c [Zx] — log(2"x) — log(Zx)© (solar abundances 
are taken from Holwcgcr 2001); 
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Table 3 

Physical parameters of the Q0347-3819 Ly« systems with metallicity gradient derived by 

THE MCI PROCEDURE 



Parameter 01 


z ab3 = 2.84829 


z ab3 = 2.89922 


(i) 


(2) 


(3) 


U 


0.1 


2.2E-2 


JVh, cm 2 


2.3E18 


1.9E18 


(7 v , km s" 1 


34.9 


54.3 


a l 


0.9 


1.3 


Z b c 


9.1E-5 


9.6E-5 


Zn 


< 2.5E-5 




Zo 


3.5E-4 




Z Si 


<2.5E-5 


2.4E-5 


[z c ] c 


-0.67 


-0.60 


[Zn] 


<> -0.5 




[Zo] 


-0.21 




[Zsi] 


< -0.14 


-0.15 


JV(Hl), cm" 2 


1.8E14 


4.1E15 


JV(Cn) 




1.1E13 


JV(Sin) 




3.6E12 


JV(Cm) 


2.7E13 


1.4E14 


iV(Sim) 




2.3E13 


JV(Civ) 


5.8E13 


3.0E13 


AT(Siiv) 


<1.3E12 


1.3E13 


AT(Nv) 


< 8.3E12 




AT(Ovi) 


5.3E13 




no. cm" 3 


2.5E-4 


2.0E-3 


(T>, K 


2.3E4 


1.2E4 


T min 


1.9E4 


9.9E3 


T max 


2.9E4 


1.4E4 


L, kpc 


3.0 


0.33 



Note. — a The internal errors of Uq, JVh, <t v , cr y and 
Z a arc ~ 15 — 20% , the column density errors arc 
< 10%, whereas no and L are known with ~ 50% accu- 
racy; b Z x = X/H; C [Z X ] =log(Z x )-log(Zx) (solar 
abundances are taken from Holwcger 2001); 
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Table 4 

Parameters of absorbers which may give rise to QSO non-DLA systems 





Type of Absorber 




Scale Size, 
L, kpc 


Mctallicity, 

X = Z/Z Q 




Mctallicity 
pattern a 




Rcf. 


1. 


HVC-typc absorbers inside 








[C/Fc] 


~ 0.2, [N/Fc] : 


- 0.6 


1,2 




galactic halos 




L < 1 


0.1 & X & 1 


[O/Fc] 


~ 0.7, [Mg/Fe]: 


- 0.2 














[Al/Fe] 


~ 0.3, [Si/Fe] : 


~ 0.3 




2. 


Gas in dwarf galaxies of 


















(i) low-mctallicity 




1 & £ & 10 


0.03 £ X £ 0.06 


[C/Fe]: 


- -0.2, [O/Fc]: 


- 0.3 


3,4 












[N/Fe]: 


- -0.2, [Si/Fe]: 


~ 0.1 






(ii) high-mctallicity 




1 < L < 10 


0.06 < X < 0.3 


[C/Fe]: 


- 0.2, [O/Fej: 


- 0.4 


3,4 












[N/Fe]: 


- -0.1, [Si/Fe]: 


~ 0.3 




3. 


Gaseous envelopes^ 


















of galaxies seen in 


















(i) Hi Lya 


10ft 


-'<L< 160ft- 1 










5,10 




(ii) Mgn AA2796, 2803 


15ft 


_1 £ L < 75ft" 1 


X < 1 








6,7 




(Hi) Civ AA1548, 1550 


100ft 


' 1 < 180ft- 1 


X < 1 








8,9 


4. 


Large-scale filaments 


















around galaxies 




L > 150ft- 1 


X < 1 


[Si/C]: 


± 0.4, [N/C]~ - 


-0.7 


10 



Note. — a The patterns for the absorbers of type 1 and 2 arc taken from Savage & Scmbach (1996) and Izotov & 
Thuan (1999), respectively, whereas the pattern for the case 4 is a characteristic chemical abundance of low-metallicity 
intcrgalactic systems (Songaila 1998). ^Following Chen et al. (1998), by 'gaseous envelopes' we mean a gaseous 
structure of large covering factor but unspecified geometry or filling factor. 
1 Savage & Scmbach 1996. 
2 Wakkcr 2001. 
3 Izotov & Thuan 1999. 
4 Mallen-Ornelas et al. 1999. 
5 Chen et al. 1998. 
6 Bcrgcron & Boisse 1991. 
7 Stcidcl et al. 1997. 
8 Chen et al. 2001a. 
9 Chen et al. 2001b. 
10 Penton, Stockc, & Shull 2002. 
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